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Abstract: In this paper the whole family of fractional Brownian motions is 
constructed as a single Gaussian field indexed by time and the Hurst index 
simultaneously. The field has a simple covariance structure and it is related 
to two generalizations of fractional Brownian motion known as multifractional 
Brownian motions. A mistake common to the existing literature regarding 
multifractional Brownian motions is pointed out and corrected. The Gaussian 
field, due to inherited "duality" , reveals a new way of constructing martingales 
associated with the odd and even part of a fractional Brownian motion and 
therefore of the fractional Brownian motion. The existence of those martingales 
and their stochastic representations is the first step to the study of natural 
wavelet expansions associated to those processes in the spirit of our earlier 
work on a construction of natural wavelets associated to Gaussian-Markov 
processes. 



1. Introduction 



The basic quantities describing movements of a viscous fluid are related through 
the Reynolds number. When the Reynolds number exceeds a critical value, which 
depends on the fluid, the average velocity of the region and its geometry, the flow 
becomes unstable and random. In his work on the theory of turbulent flow Kol- 
mogorov proposed a model which assumes that the kinetic energy in large scale 
motions of a turbulent flow is transferred to smaller scale turbulent motions. At 
smaller scales the Reynolds number associated with that region is reduced. When 
the Reynolds number of a region falls below the critical value of the region, turbulent 
motion stops and the remaining kinetic energy is dissipated as heat. Kolmogorov 
assumed that smaller scale turbulent motions can be described by a random field. 
Kolmogorov modeling assumptions produce the random field Bh that is self-similar 
with stationary increments and the second moment of its increments is of the form 
E \Bu(t + s) — Bh (s)| 2 = c|i| 2ff . Analyzing the properties of those fields [1,0] 
Kolmogorov obtained a spectral representation for the fields with stationary incre- 
ments. 

In 1968 Mandelbrot and Van Ness interpreted the nonanticipating representation 
of B H , 

B H (t) - B H (s) = C J (ft- x) + ~ i - (s - db(x), 
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with respect to an orthogonal white noise db, obtained earlier by Pinsker and Ya- 
glom, as a fractional integral, and called Bh the fractional Brownian motion in 
the case when b is a Brownian motion. Most of the above information is quoted 
from Molchan [3| • His paper is an excellent short review of the history of fractional 
Brownian motion before 1970. The index H is sometimes called the Hurst exponent, 
after the British hydrologist H.E. Hurst, who studied the annual flows of the Nile. 

Fractional Brownian motion has applications in financial mathematics, telecom- 
munication networks, hydrology, physical and biological sciences, just to mention a 
few. 

An appropriate representation of fractional Brownian motion is important for 
analyzing its properties and for computational purposes. Meyer, Sellan and Taqqu 
developed a method for constructing wavelet expansions of fractional Brownian mo- 
tion [12| • Their construction was based on fractional integration and differentiation 
of wavelet expansions of a Gaussian white noise and encompasses a large class of 
wavelets. 

We have found an iterative method for obtaining an orthogonal expansion of 
a Gaussian Markov process [B| directly from its covariance function. It turns out 
that our method produces a wavelet expansion if time is measured by the natural 
measure associated to the process. It is therefore natural to ask if it is possible 
to construct a "natural" wavelet expansion associated with a fractional Brownian 
motion (fBm in short) in the spirit of our work in 5]. The main properties used 
in our construction of wavelets associated with the Gaussian Markov processes 
were the invariance of the processes on pinning (stays Markov), independence of 
associated pinned processes, and the existence of associated martingales. 

The first step toward finding natural wavelets for fBm is to investigate invariances 
for the entire class of fractional Brownian motions. That means considering the 
family of processes (Bh) 0<h<1 where Bh = (Bh (t)) teK is a fBm. Let us point 
out that it is not straight forward to prove that, for the fixed Hurst index H, the 
covariance function, 

(1) EB H (t) B H (s) = \{\t\ 2H + \s\ 2H \t - s\ 2H ), 

of fBm is indeed a positive definite form . It uses tricks and the fact that the char- 
acteristic function ip of an 2H stable random variable is given by ip(t) — e~ a l*l for 
some constant a. In the search for invariances relevant to achieving our goal, that 
is to obtaining natural wavelet expansions for fBm, we have constructed a Gaus- 
sian field, indexed by K x (0, 1) , that encompasses at once all fractional Brownian 
motions for all H. More precisely, this field has the property that, when the second 
parameter, H , is fixed, the resulting process is an _ff-fBm. We started our construc- 
tion by using already mentioned nonanticipating representation of fBm, which, for 
< H < I, can be written as 

(2) B H {t) = —j- n / (t-x) + -(-£c) + 2 dW(x) 

for t e R, where (W (t)) teR is a Brownian motion. When we started the compu- 
tations our hope was that the covariance of our Gaussian field would resemble the 
covariance of fractional Brownian motion (equation {!])), with 2H substituted by 
H + H modulo a function of H and H . This does not turn out to be the case. The 
resulting covariance was more complex. Its form revealed that the nonanticipating 
representation may not be the most " efficient" representation of fBm. Section [3] of 
this paper contains all computations relevant to finding that covariance. 
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Although the Gaussian field developed in Section[3]did not meet our objective, its 
structure revealed what intrinsic dependences of the nonanticipating representation 
of fBm are, and what should be done to construct a new field having, for our 
purposes, the "right" covariance, namely the covariancc 

(3) EB H (t) B H , (s) = a H H , (\t\ H+H ' + \s\ H+H ' \t s\ H+H ' ). 

We have succeeded to obtain such a field. Furthermore we have found a nonan- 
ticipating representation of {Bh (i)lteR i?e(o 1) ' wmc h, when H is fixed, equals in 
distribution to the standard nonanticipating representation of fBm. Those results 
are presented in Section 21 

Our new field has the property that when H + H =1 the right hand side of 
© becomes the covariance of a Brownian motion. In that case we call Bh, B h i a 
dual pair. A particular property of a dual pair is that it generates two martingales, 
one driving Bh and the other B H > . In [J] we have obtained stochastic integral 
representations for those martingales. It turns out that our representation coincides 
with the fundamental martingale of a fBm discovered by Molchan and Golosov 
dUU^]. Our work on that subject is in Section [5l 

Levy Vehel and Peltier [23[ and Benassi et al. [IJ have independently introduced 
a multifractional Brownian motion. In Section[6]we clarify the relationship between 
multifractional Brownian motion and the fractional Brownian fields introduced in 
Section [3] and [4] Furthermore we point out and correct an error in the covariance 
of the multifractional Brownian motion obtained from the nonanticipating moving 
average of fBm and show that in fact the processes of Levy Vehel and Peltier [23j 
and Benassi et al. 0] are not the same, as it has been claimed in Cohen 

For readers convenience some basic facts about fractional Brownian motions and 
the notation that is used in the paper is included in Section [2] 



2. Notation and preliminaries 

This section is a brief overview of some properties of the fractional Brownian motion. 



The material presented here can be found in Taqqu [22j, Molchan [14j. chapter 7 
in Samorodnitsky and Taqqu 21 1 and Dzhaparidze and Van Zanten [6|. All the 



processes and fields considered will be real valued. 

The following proposition summarizes some properties of ii-self-similar processes 
with stationary increments and finite second moments. 

Proposition 1. If Z = {Z t } tgR is H -self-similar with stationary increments and 
finite second moments then 

1. Zq = with probability one. 

2. If H ^ 1, then EZ t = for all t e M. 

3. H < 1. 

4. The covariance of Z is given by 



(4) E (Z s Z t ) = [\s\ 2H + \t\ 2H -\s- t\ 2B 

As already mentioned in the introduction Kolmogorov described an _ff-self- 
similar process with stationary increments and its covariance (H in 0]. A centered 
Gaussian process {Z t } tem with covariance ((4|) is called a fractional Brownian mo- 
tion (fBm) with index H. If E (Zf) = 1, the process is a standard fBm. Since 



so 
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Bi (t) = 



the distribution of a Gaussian process is entirely determined by its covariance, ([4]) 
implies that fBm is iiZ-self-similar and has stationary increments. Conversely, an 
ii-self-similar Gaussian process with stationary increments is a fBm. For H = i 
a standard fractional Brownian motion is a standard Brownian motion on the real 
line, that is, 

B 1 (-t) t < 
B 2 (t) t > ' 

where {i? t 1 } t>Q , {i? t 2 } t>Q are two independent standard Brownian motions. 

Unless otherwise stated all Brownian motions and fractional Brownian motions 
appearing on this work will be assumed to be standard. 

Using Kolmogorov-Centsov criterion for Holder continuity of a process, for every 
7 in (0, H) , an i/-fBm has a continuous modification that is locally Holder contin- 
uous with exponent 7. Throughout the paper we will assume that fBm is such a 
modification. 

Let {Bt} t(£R and {Wt} t€R be Brownian motions on the real line. Then the process 
defined by the moving average stochastic integral 

1 f 00 1 1 

(5) B H (t) = — / (t - x)"~* - (-x)"' 1 dB x for t e K 



ch 

is a fBm, where 



(6) c H 



V / T(2H +l)sin {ttH) 



and x" = if x < and x°L — x a if x > 0. This representation is called nonantic- 
ipating, since for t > 0, it involves only integration over (— 00, i]. Another moving 
average representation of fBm for H ^ ^ is given by 

(7) W H {t) = -j- [ \t-x\ H ~^ -\x\ H ~^dW x forieIR 

">H J -co 

where 

(8) d H = 




v / T(2H+l)]j V sin(7ri?) 
and when H = | by 

(9) Wi (t) = J- / log f r— !— r") - log ( dIU x for * S 



di ./_^, V k — x 



where di = ir. This representation is called well-balanced. 

The odd and even part of fBm play an important role in this paper. If {Zh (^)} te R 
is a H-iBm, then the odd and the even process of fBm are defined by 

Z° H (*) = \ (Zh (t) - Z H (-*)) , and (t) = ~ (Z H (i) + Z ff (-*)) , tgK 

respectively. It is straightforward to check that the processes Z° H and Z H are inde- 
pendent and that their covariances are given by 

(10) EZ° H (s)Z H (t) = ±(\s + t\ 2H -\s-t\ 2H 
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and 

( u) EZ . M ^ w= kr±g_ k±ir±^f! 

for s, t in E. Both ZJ^ and are TJ-self-similar processes and by assumption have 
continuous sample paths with probability 1. Clearly, each path of the process Z° H 
(Zfj) is an odd (even) function on a set of probability 1, and so it suffices to consider 
these processes for t > 0. For H = | the odd and even part of fBm are Brownian 
motions (up to a constant multiple). 

Moving average representations for the odd and even part of fBm can be found 
in Dzhaparidze and Van Zanten [6j], Nuzman and Poor [17| . 



3. Dependent fractional Brownian field 

In this section we will assume that each element of the family {Bh (t)}t^s. He(a l) 
is represented by a nonanticipating moving average stochastic integral ([5]) respect 
to the same {B t } tm . For < H, H' < 1 set 



_ y/f (2H + 1) sin (ttH)^T (2H' + 1) sin (irH') 
c h,h' — Z ■ 

Theorem 2. Let K be the covariance EB H (t) B H > (s) . If H+ H £ 1 then 



A =c HH ,r(- [H + H'yj jcos((V -if) |) cos^H + H'^ " 



x [\t-s\ H+H -\t\ H+H -\s\ H+H 
- sin ( [h' - H^j |) sin (Jh + II 



')!) 

x ( sgn (t) \t\ H+H ' - sgn (s) \s\ H+H ' - sgn (t -s)\t- s\ H+H ' 



If H +H' = 1 andt^ s, t^ 0, s ^ 0, then 

K = c HM , {cos ((V - H) 1 ) | (\t\ + \s\ - \t - S |) 



(13) -sin((V-i/)^ 
x (t log |t| - slog \s\ - (t-s)log|f-«|)}. 

If H +H' = 1 andt = s, t ^ 0, i/ien 

(14) X = c HjH , cos ((if' - F) I) | {|t| + | S |} . 

If H +H' =1 andt = or s = 0, then 

(15) K = 0. 

Fourier transform is the main tool in proving the theorem. Before proving it we 
will establish some technical results. 

Let / (£) = e %( - x f (x) dx be the Fourier transform of a function that belong 



-oo 

- 2 



to L 1 (E) n L 2 (E) With this convention v2tt ||/|| 



for all / £ I 2 (I 

2 V 
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Lemma 3. For < H < 1 and t <G K 



w/iere 

(18) «n ^ - J 11 



> 



Proof. Equation (|16p follows from the theory of fractional integration and differen- 
tiation, see for example Lemma 1 and Lemma 3 in (20j . Let 

r(i/ + i) 

then the identity 

rWTf) = f -' ( - x) 

is a consequence of the relation x°t. = (—x)_. Equation (jTTJ) now follows from 

dm). □ 



Lemma 4. For s,teK and < H,H < I let 



, . roo sin 2 ( | ) + sin 2 ( f ) - sin 2 ( ^ 







IfH + H ^ 1 then 



(20) 



H + H ) ) cos ( ( H + H )f 



2 



and if H + H =1 i/ien 

j 1 = I(|t| + H-|t- a |). 

Proof. The proof is trivial if either s = or t = 0, so assume that s ^ 0, t ^ 0. 
Formula 3.823 in Gradshteyn and Ryzhik [7[ states 

r°° r fu) cos (^) 

(21) y a;^ 1 sin 2 (ax) dx = 2A1+1 \ 2 for a > and - 2 < Re (p) < 0. 
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When H + H' ^ 1 set fj, = - (^H + H'^J . Observe that sin 2 (x) = sin 2 (|a;|) . Apply- 
ing the identity (f2Tj) to the right hand side of (fT9|) when s = t gives 



h = -r (-H - if') cos (Jh + H >S j ~) \t 

and when s^t 



-H-H 



cos((ff + ff') |) (3(H + H) 



x{\t-s\ H+H -\t\ H+H -\s\ H+H 



In the case of H + i? = 1 we will use formula 3.821 (9) from Gradshteyn and 
Ryzhik 0, 



sin (ax) an 

dx — — for a > 0. 



to conclude that when s ^ t then ii = 7r/4(|t| + |s| — \t — s\) , and when s — 
t then Ii = 1/2 1*| vr. □ 

Lemma 5. For s, t S M and < H, H' < 1 Ze£ 

de/ r sin (g (t - 3 )) + sin (gQ ~ sin (tQ 
[ ' 2 Jo £1+H+H' 

If s = t or s = or t = £/ien /2 = 0. Otherwise, if H + H ^ 1 f/ien 
(23) 



7 2 = r ( - ( .ff + if ) ) sin ( ( H + H ) ~ 



2 

x |s 5 n (i) - sgn (s) \s\ H+H ' - sgn (t - s) \t - s\ H+H ' 

and if H + H =1 then 

(24) I 2 =t log \t\ - slog \s\ - (t-s)log|i-s|. 

Proof. The proof is trivial if either s = 0orf = 0ors = i, so assume that s ^ 0, 
t 7^ and s 7^ i. Formula 3.761 (4) in Gradshteyn and Ryzhik Q reads: 



(25) 



/ x^ 1 sin (ax) dx = sin (^-) for a > and < |Re (/i)| < 1. 

Jo o/ 1 V 2 / 



Set p = — (H + Hj , and observe that sin (£x) =sgn(a;) sin (£ Applying (f2"5")l 
to the right hand side of (|22[) in the case when < H + H < 1 yields (|23[) . In the 
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case when 1 < H + H < 2, by the dominated convergence theorem, it follows that 
= f b sgn (t - s) sin (g |f - sj) + sgn (a) sin (|s| Q - sgn (t) sinjt] £|) 

b— *oo 

l /" b 

= lim + wTTTr / ( s # n " s ) I* " s l cos (£ F ~ S D 

a^0+ H + H J a 
b — >oo 

+sgn (a) \s\ cos £) - sgn (i) |t| cos (\t\ £)) f ff ^' ^ 
= (H + H')(H + H'-l) l (-s 9 n(t-s)\t-s\ 2 ^\t-s\) 
-sgn (s) |s| 2 sin + .sgn (t) |t| 2 sin (|t| £)) f 1 "*-^ d£, 

where the last two equalities are result of integration by parts and the fact that 
the boundary terms converge to zero as b — > oo, and a — * + . Applying (|25[) with 

fj, = 2 — ^Jf + i?'^ and using aT (x) = r (x + 1) the equation now follows 

readily. 

Let us turn to the case H + H =1. For x £ [1, 1.5] set 



n 



,. ( r) = i sgn {t - s) sin (gjf--a[) + sgn (s) sin (|s| Q - sgn (t) sin g) ^ 



When a; = H+H the function / equals to the right hand side of (|23l) . The integrand 
that defines / is bounded by 



\sgn(t-s) sin(£\t-s\)+sgn(s) sin(|s|g)-sgra(t) sin(| | £ £ (Q j\ 



which is an integrable function. By the dominated convergence theorem / is con- 
tinuous on [1, 1.5]. Finally lim^i / (x) establishes (f24|). where we used (|23|) and the 
gamma function property yT (y) = T (y + 1) to rewrite / (x) for x S (1, 1.5] as 

^ = (-i) 2 (i -x) sin (t) ^ S5n ^ ~ S5n ^ ' S ' :E ~ S5n ^ - s ) ^ - s ri 

and L 'Hospital rule to compute the limit 

lim Sgn ® ~ Sgn ^ ' S ' X ~ Sgn ^ ~ ^ '* ~ 



sj.1 (1 - x) 

We have prepared the groundwork to prove Theorem [2l 
Proof. By the Ito isometry 



□ 



K = 

ch 



l — jT ((t - a;)f-* - (-a)?" 4 ) ((« - *)+ "* - (-*)+"*) 
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and by Planchcrel identity and Lemma [3] 



r(jj + i)r(V + i) 



a- = — * '- \ —— {ity^-^ —— -^de 

2ttc h c h > J_ 0o i£ i£, 



io |£| J 

Substituting £ = — £ in the first integral of the last equality above and then com- 
bining the two integrals yields 

r (g + |)r( g ' + |) / f~ ( e «e _ i) ( e -w _ i) \ 

Using Euler's formula on e l ( H H )~ and (e 1 *^ — l) (e _ls, » — l) and the identity 
sin 2 x = 1-c ° s 2x we obtain 

Re f e <-'-)f r o^-iH^-iU 



7T 



DC- 



sin 2 (f)+sm 2 (#)-sin 2 (^) 



(27) =cos^ -H)-)2j o — ^ > ^ 

- sin ((H - H) I) r toMt-s))+J*W-**(X)^ 
Jo C 

Observing that 

r(H + h)r(H' + l 



(28) 



yT (2H + 1) sin (nil) y/T (2H' + 1) sin (ttH') 



chc h '1t 



the expressions for E£># (t) B H > (s) now follow from equations (f2"o| , (|27p , (|28p and 
Lemmas HI [5] □ 

We will call a centered Gaussian field {-Bj/ WltgR He(o i) w ith the covariance 
given by Theorem [2] a dependent fractional Brownian field and refer to it as dfBf. 
The rest of the section elaborates on a property of the field that justifies that name. 

Let {B° H (t)} te 

[o,oo),He(o,i) anc ^ {^h (*)}te[o,oo),He(o,x)' ^ e tne oc ^ an< ^ even P art 
of the dfBf {B H {t)} teR He(0A) , that is 

B° H (I) = and B% (t) = g^l + g^M, , > . 

For + H' ^ 1 set 

y/T (2H + 1) sin (ttH) y/T (2H' + 1) sin (ttH' ) 



a H,H' ~ ^ 

(29) 



7T 

7T\ //„ '\ 7T N 



x T ( - | II - II ) ) cos ( ( // -H)-)oob([H + H , } 
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and for H + H' = 1 



(30) a H,H' = VriZH + 1)T(3- 2H) sin 2 (irH) . 

Theorem 6. The covariance of the odd part and the even part of dfBf are 



(31) EB° H (t) B° H , (s) = a HtH , ^ 



and 

EB e H (t)B e H , (s) 



(32) 



\t\ H+H +\s\ H+H \t-s\ H+H +\t + s\ H+H 



l H,H 



respectively. 

Proof. The covariance of the dfBf (Theorem \2§ , is of the form 



EB H (t)B H , (s)=f(H,H') (\t\ 

(33) / V 

+ g(s,t > H,H 



H+H +{sl H + H _ ]t _ s] H+H 



It is a matter of straightforward computation to check that in both cases 

EB H (t) B H , (s) + EB H (-t) B H > (s) and EB H (t) B H > (s) + EB H (-t) B H > (s) 

the g function cancels. The result of the theorem follows by simple algebraic ma- 
nipulation of the first part of the right hand side of (1351) only. □ 

The Gaussian fields {B° H (t)} te[0 oo) He(01) and {B C H (t)} te[0 oo) He(01) with co- 
variances given by Theorem [SJ will be called the odd and the even fractional Brow- 
nian field respectively. It is very simple to check that for every a > 

{B°h ( a ^}te[o,oo),_f/e(o,i) /= i -" B °h ( a< )}t S [o,oo),ffe(o,i) 

and 

{B e H (ot)} tE [o )00 ),He(o,i) = { aIiB H ( at )} t e[o,oc),He(Q,i) ' 

where ^ = d ' indicates the equality of finite dimensional distributions. 

Given a fBm its odd and even part are independent processes (indexed by 
t). However, this is not the case with the dfBf {Bh {t)} teW Hel ^ The fields 

{B° H (*)}te[o,oo),ire(o,i) and i B H (*)}te[Q,«0,ffe(o,i) are not independent. For exam- 
ple if H + H ' = 1 then 



EB e H (t) B° H , (s) = -y/T{2H + l)sin( 7 r J ff) 



x \/T(2H' +l)sin(7rif') si n((-H"' - H J | 

(f + «)log|f + «|-(f-«)log|f-«| 



x i slog |s| - 
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which is clearly not equal to 0. That is the reason for calling that field the dependent 
fractional Brownian field. 

Another glance at the computation of the covariance EB' l H (t) BP , (s) , i,j S 
{o, e}, reveals that when i = j the g part of (|33[) cancels out while in the case when 
i ^ j the first part of (f33|) cancels out leaving the g part. Therefore the existence 
the g part in (|33j) is the reason for dependence between {-B|j (i)}*eR _hg(o i) an< ^ 
{£?f^ (i)} t6K // S (o i) ■ S° it is natural to search for a method of creating a fractional 
Brownian field that would be of the form ([33]) with 5 = 0. One way of attacking 
that problem is a direct verification of positive dcfiniteness of such a form, a very 
unattractive task. In the next section we present a straightforward construction of 
the field with the desired covariance. 



4. Fractional Brownian field 

The last remark of the previous points the direction for constructing a fractional 
Brownian field {Bh (OlteR He(o l) with the covariance of the form of (|33| with 
g = 0. The new Gaussian field contains all fractional Brownian motions too. 

Theorem 7. Let B = {B H {t)} tmMe {o.i) and W = {W H (*)} ieR ,H e ( ,i) be two 
dfBf generated by two independent Brownian motions {B t } teR and {Wt} t€R respec- 
tively. Let {B H (i)} te [ oo) He(o l) >* = ^ e ^ e oc ^ an< ^ i = e be the even part of 
B, and let {W H (i)} tg j ^ hg(o i) ' ' = ^ e ^ e an< ^ * = e ^ e even P ar ^ °f W ■ 
Then the fractional Brownian field {Zh (i)} tgK ue(0 l) defined by 

B e H (t) + W H (t) fort>0 
B e H (~t) - W° H {-t) fort<0 



\t\ H+H ' + \s\ H+H ' - \t-s\ H+H ' 



(34) Z H (t) : 
has the covariance 

(35) EZ H (t) Z H , (s) = a H M 



where an.H' * s given by equations (|29[) and (|30[) . 

Proof. The proof follows from (|32|) . ([31]) and independence of {-Bt} tgR and {Wt} tgR . 

□ 

We will call the process {Zh (t)} teR ne(o l) fractional Brownian field (fBf in 
short). Note that for any tel, 

Z H (t) = ~ (B H (t) + W H (t) + S H (-t) - W H (-t)) , 

and that an d J B B 1A-w H {t) \ are twQ independent 

I v 2 J teR, L v 2 J tea,j?e(o,i) 
dfBf. Consequently 

(36) Z H{ t) = ^miEUl 

v2 

where {Xp (t)} tgR ne(o l) anc ^ WlteR ffe(o i) are ^ w0 independent fractional 
Brownian fields, that is {Zh (i)ltgR He(Q l) * s a P ro P er b/ symmetrized dfBf. 
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Proposition 8. Let {Bt} teR and {Wt} t£R be two independent Brownian motion 
processes on the real line and let {Zh {t)}t^R he(o i) be defined by (|34|) . Then 

Z H (t) = -=— / (t - xf + -> - dB x 
I f°° rr i ft i 

Proof. Follows directly from (f3"6")) and ([5]). □ 

A fBf has the same self-similarity property in the time variable as the odd and 
even fractional Browian field, namely for a > 

{Z H (at)} tm He{a l) = ' {a H Z H (at)} teR He{01) . 

Moreover, the stationary in the time variable of increments of the fBf easily follows 
from Theorem [3 that is 

{Wh (*)} teM r e (o,i) f '=' {Wh (t + S)- W H (*)} teB)J?6(0)1) 

for any S. 

An immediate consequence of the covariance structure of a fractional Brownian 
field is that when H + H' = 1 then 

,37) E^W^WJ-^IW^^i'-^-. 

That property leads to a construction of martingales associated to fractional Brown- 
ian motions. The methodology of the construction is the subject of the next section. 



5. Duality and fundamental martingales 

In what follows it is assumed that {Bu (£)} tgK hg(o i) ^ s an ®^ Whenever H+H' = 
1 we will call Bh , (or B H ,B° W or _B|^, -Bfp) a dual pair. Dual pairs have unique 
properties. They generate martingales associated in a natural way to fractional 
Brownian motions Bh and B^* . The construction and explanation of the nature of 
those martingales is the subject of this section. 

Every fBf is a sum of an even and an odd part of two independent dfBf 's. For that 
reason it suffices to construct martingales, and M^, adapted to the filtrations 
of the odd and even part of {Bh (£)} (6 r respectively. The filtration generated by 
Mfj (Mfj) coincides with the filtration of the odd (eve n) p art of fBm, and for that 
reason, following the terminology used in Norros et al. [161 ] to describe a martingale 
for the fBm originally discovered by Molchan and Golosov (see Molchan 14|), we call 
Mjj (Mjj) a fundamental martingale for the odd (even) part of fBm. Furthermore 
we derive a stochastic integral representation for those martingales. In a similar 
fashion this was done in Pipiras and Taqqu [TU and Pipiras and Taqqu [3] for the 
fractional Brownian motion. 

For i E {o, e} set 

Tf* = a (B l H (s) : < s < t) and G? = span (B l H (s) : < s < t) , 
where {B° H (i)} t>0 and {B e H (t)} t>0 are the odd and even part of {B H (t)} teR . 
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For t > 0, i e {o, e} and H + H' = 1 define 
(38) Mh (t) = E (B l H , (t) | F?*) . 

Theorem 9. {Mfj (t)} t>0 and {Mfj (t)} t>0 are H- self- similar Gaussian martin- 
gales adapted to the filtration {ff'°}oo and {F t e } t>orespectively . 

Proof. It is enough to verify the statement for M\\ only, because the verification 
for Mfj is similar. By construction Mfj is a Gaussian process. It follows from (|3Tj) 
that for s <t, 

B° H , (t) B° H , (s) ± Gf °, 

which implies 

E (Mjj (t) | = E (l (fl^ (t) | *f | = E (fl^ (t) | ■ °) 

= E (B° H , (t) - (s) | .Ff >°) + E (B° H , (s) | 
= + M^(s)=M^( S ). 

By iJ-self-similarity property of the odd part {B° H (t)} te r ^ ffG ( l) °^ ^ ne ®^ 
{^h (*)} tg ^H 6 (o,i) the field ^ (*)}ie[o,oo),He(o,i) defined by 

Z£ (i) = a"^ (of) 

is an odd fBf and, therefore for H + if = 1, 

{E (2^ (t) | a (Z& (r) : < r < i)) } t > ^ {M° H (i)} t > . 
Furthermore 

(j (Z° H (r) : < r < t) = cr (a~ H B° H (ar) : < r < f) = cr (ar) : < r < t) 
= o-(B° H ( S ):0< s <at)=^°. 

Therefore 

E {Z° H , (t) \ a(Z° H (r):0<r<t))=E (a~ H B%, (at) \ , 

which concludes the proof. □ 

So far we have shown that M\\ and Mfj are 7J-self-similar Gaussian martingales. 
By construction, for i £ {o, e}, M l H (t) is an element of G^' 1 , and therefore it may 
be possible to express it as a stochastic integral, up to time t, of B l H . In the case 
H = I this is trivial, since then i? = | and therefore -B^ = -B^., is a constant 
multiple of Brownian motion and M l H (t) — B\j (t) . The case H ^ h has been 
solved in our paper 0]. We state the result below without proof. The supporting 
materials are too long for the present paper. It should also be mentioned that the 
natural filtration of the martingale Mjj coincides with the natural filtration of the 
process B % H 0]. 



Theorem 10. Let H e (0, 1) \ Ift > then 

M° H (t) = Y ^^j\t>- S >)i- H dB° H (s) 
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M% (t) = 



r(l-H)J ds 



(x 2 - s 2 ) 2 " dx ) dB% (s) , 



where 



a H = 



2 2g -Vr(3-2iJ) sin(Trg) 
T (I~ H ) y/T(2H + l) 



6. Remarks on multifractional Brownian motions 

Levy Vehel and Peltier [23| and Benassi et al. [2| have introduced independently, 
multifractional Brownian motion. In this section we will clarify the relationship 
between multifractional Brownian motion and the fractional Brownian fields intro- 
duced in Sections [3] and 21 Additionally we will point out an error in the covariance 
of multifractional Brownian motion obtained from the nonanticipating moving av- 
erage representation of fBm which shows that in fact the processes of Levy Vehel 
[23| and Peltier and Benassi et al. Q are not the same, as it has been claimed in 
Cohen Q. 

Let {Ws} tgK be a Brownian motion. For t > Levy Vehel and Peltier [23| called 
(39) X t = j— i ■ n / Ut-s)^ - (-a)?*-*) dW Sl 



(H(t) 



2/ 



where H : [0, 00) — > (0, 1) is a deterministic Holder function with exponent f3 > 0, a 
multifractional Brownian motion. This process is introduced as a generalization of 
fBm that has different regularity at each t, more precisely, if < H (t) < min (1, (3) 
then at each to the multifractional Brownian motion has Holder exponent H (to) 
with probability 1. It is clear that if H (t) = H for some < H < 1, then {X t } t>0 
is a (nonstandard) H-fBm. 

Benassi et al. [2| have introduced the process 

(40) Yt = Lwm M " 

where in "some sense" the random measure dW is the Fourier transform of dW 
and for g,h 6 L 2 (R) it satisfies 



El/ g (0 dWs ^ h (0 dWt j = f g(0h 

(see section 7.2.2 in Samorodnitsky and Taqqu If H (t) = H, for some 

< H < 1, the process {V t } 4>0 is an H-fBm, because the right-hand-side of (jlCJ)) 
reduces to the well-known harmonizable representation of fBm. The result concern- 
ing the Holder exponent for {X t } t>0 holds for the process {V t } t>0 too. Although 
the process {5 / t} t >o is called multifractional Brownian motion, we will refer to it 
as a harmonizable multifractional Brownian motion to emphasize the differences 
between the two processes. 

In Q Cohen states that both multifractional Brownian motions {V t } t>0 and 
{Y t } t>0 if normalized appropriately are versions of the same process. More precisely 
the following is stated in [3] (as Theorem 1): 
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The harmonizable representation of the multifractional Brownian motion: 
(41) / -TTW^dW,, 



is almost surely equal up to a multiplicative deterministic function to the well bal- 
anced moving average 



(i 



t _ s \H(t)-h _ \ s \ H m-?) aw. 



When H (i) = i, — s\ H ^ 2 — |s|^*' 1 2 ^ is ambiguous, hence the conventional 
meaning 



^ s ' -' s '° = lo HF^J- log VN 

is to be used. Conversely, one can show that the non anticipating moving average 
(42) f ((t-sf + {t) -* -{s)l {t) -*)dW s 

is equal up to a multiplicative deterministic function to the harmonizable represen- 
tation 

r ui _ i — 
-dW ( . 



1 \ , / 1 



Kl £|£| 



H(ty 



Hence the mfBm given by the non anticipating moving average (|42|) has the same law 
as the mfBm given by the harmonizable representation (|4ip up to a multiplicative 
deterministic function. 

The arguments used in the proof of the of Theorem 1 in Cohen Q are based on 
the fact that the Fourier transform of x -> \t- x\ H{t) ^^ - \x\ H{t) ~^ for H (t) ^ 
\ and of log(-u^r) — log(A) for H (t) = | equal, up to a multiplicative con- 

stant, to £ — > e . , — =— ; and an incorrect statement that the Fourier transform 

|f|2+ H <*> 

1. Hit) — 

of x — > (t — ie) , 2 — (—a;) i 2 equals, up to a multiplicative constant, to 
£ — > — e t . The equations (116D and ([TBI) are the correct expression for that 

Fourier transform. Consequently the last two statements of the above Theorem 1 
in Cohen are incorrect. In order to see why, consider two multifractional Brownian 
motions 

i r 00 i i 

X t = — (t- xf + {t) ^ {-xf + {t) ^ dB x 

c H(t) J-oo 

and 

y f l« - - ^ fOT ' * W * 3 

* 1 ^/^log^-log^)^ *(*) = §' 

where £ > 0, Ctf( t ) and d H ^ arc defined by ©, © respectively, di — it, and 
{^*iteR ano - {^iteW-i 9,16 Brownian motions. According to the last statement of 
Theorem 1 in Cohen |3[ there is a deterministic function f t such that the processes 
{X t } t>0 and {ftYt} t>0 have the same law. The chosen normalization assures that 
E (Xfj = E (Y t 2 ) for'all t > 0, implying that |/ t | = 1 for all t such that t > 0. It 
follows that |E(X(X S )| = |E(y t y" s )| for all s,i. It is clear that the process {X t } t>0 



92 



V. Dobric and F. M. Ojeda 



can be obtained from a dependent fractional Brownian field {Bh (t)}t>o ffe(o 1) as 
{Bn(t) {t)}t>o- Similarly, the Gaussian field {Wh (t)} t>Q H e(o n denned for H ^ \ 
by equation ([7]) and for H = \ by © gives {Y t } t>0 via {Wn(t) M} f>0 - Since the 
last statement of Theorem 1 in [3| is supposed to hold for every Holder function t — > 
H (t), that statement holds if and only if the Gaussian fields {-E>ir(i)}tg[o,oo),.f/e(o,i) 
and {Wh (i)}te[o,oo),_ffe(o,i) have the same covariance in absolute value. The covari- 
ance of {Bh {t)}te[o,oo),He(o,i) * s gi ven by Theorem [SJ Proposition [TT1 right below 
gives the covariance of {W07(i)}te[o.oo)..£re(o.i)- The proposition shows that if H ^ -| 
and H+H = 1 the covariance EW# (t) W H > (s) is a multiple of sAt. From Theorem 
[2] we see that this is not the case for Ef?# (<) B H > (s) . 

Proposition 11. The covariance of the Gaussian field {Wh (t)}t^m He{o l)> defined 
for H 7^ i by equation ([7]) and for H = ^ by © is 



I I d 2 > I ,,,H+H . I ,H+H I, ,11 + 11 

JJ ^2 H + H 



where dn is defined by (0) and 

(43) fc H = -2r \H + i ) sin [ ( H - - ] - ) for H ^- and ki 



27 VV 2) 2) J '2 
Before proving Proposition 111! a few technical results are needed. 
Lemma 12. Let f t i (x) = log j^z^y — log tn- /or x.tel. TTien 

e it? _ i 



/a (0 



Proo/. Suppose that £ 7^ 0. Since f t i e L 1 (R) n i 2 (K) the Fourier transform of 
f t 1 can be computed as 

(fl - £ (>» g ^ - * i,) e-^ . £ log 
Substituting w = x — | yields 

and after the substitution u = —v on (—00, 0) to 

, , r°° ( \ u 4. i| \ 

/ t 1 (0 = e^i ^ log f 1— -|l J sin «) du. 

Formula 4.382 (1) in Q reads: 

/ log I : I sin (bx) dx — — sin (ao) for a, 6 > 0. 

If i > set a = 5 and 6 = |£| and use that sin = son (£) sin (z |£|) for z > to 
get 

^(£) = e^ S m(§ 
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Standard trigonometric identities sin (2a) = 2 sin (a) cos (a) and sin 2 (a) = 1 c ° h ( 2a ) 
complete the proof for t > 0. The proof for t < is similar. □ 

Lemma 13. Let f t , H (x) = \t - xf~^ - \x\ H ~^ for x,t G R and H e (0, 1)\ 
Proof. Rewriting / t # as 



27 2/ |^ 



and applying equations l|ldp and 1)170 it follows that 

n erf? - 1 - -(*-*) _ r 4. ^ ----1 r_*r(*- 



(0 = r ( + ,j — - r I * + - 2 ) — HO 



and from ([T8l) that 



(^"(^-i) - He)-M) = 2 * sgn (-0 sin ((ff - 



1\ 7T 

2 



Therefore 



M _ 












sin 


:(- 











1\ /V 1\ tt\ e it( - - 1 

= 21 \H 



2) \\ 2)2) 

□ 

We are now ready to prove Proposition[TlJ The idea is to use the fact that up to a 

multiplicative constant f t ,H (£) f s h' (0 equals / H+H ' (£) / H+H > (£) and that up 

to a multiplicative constant the integral over R of the later is the covariance of an 

J 1 -fBm. This argument is used in Ayache et al. [1[ to compute the covariance of 
a multifractional Brownian motion given by (|4TJ|) . In Ayache et al. |l| it is also erro- 
neously claimed that the covariance of the multifractional Brownian motion given 
by (|39p is the same, if properly normalized, as the covariance of the harmonizable 
multifractional Brownian motion given by (|40[) . Their proof is based on the last 
statement of Theorem 1 in Cohen In section 3 of Lim and Muniandy [lH, the 
authors give another incorrect argument about the equivalence (up to a determin- 
istic multiplicative function) between the harmonizable multifractional Brownian 
motion (|40[) and the nonanticipative multifractional Brownian motion (|39p . 

Proof of Proposition By the Plancherel identity 

1 f°° 

EW H (t) W H , (s) = —— / f ttH (x) f s H , (x) dx 
(44) " ^ ~ oo 



2tt dud H ' J_ 



/GO 
-CO 
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where ft,H> when H = i, is the same as in Lemma [T2l and, when H ^ |, as in and 
Lemma [T3] From Lemma [12] and Lemma [13] it follows 



r°° , _ -^=^z f (e u £ - l) (e ls i - 1) 



Let H = H+H . Observe that H S (0, 1) and so 



r 00 - — ■ d 2 „ k 2 H f (e lte - - 1) (e ,s « - 1) 

cP 2ir f 00 

= ^H^H'^-n- \ ft,Ho ( x ) fs.Ho ( x ) dx , 

where the last equality follows from Lemma [T2l Lemma [T3] and Plancherel identity. 
Hence 

d 2 

— "' -s) 



/oo - (jf 

ft,H (0 f, tH > (0« - k„k H , -§±2ttEW Ho (t) W Ho ( 
-oo K H 



k H k H >j§ 2 -2'K 
k h 



\t\ 2H " + \s\ 2H " - It - s\ 2Ho 



Hence (|44|) becomes 

fefffc^ d 2 Ho ( \t\ 2H " + \s\ 2H ° - |t - s| 2H ° ' 



IW H (t) W H , (s) = 



d H d H > k 2 Ho y 2 J' 

This finishes the proof. □ 
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